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HiggsBounds is a computer code which tests the Higgs sectors of new models against the cur- 
rent exclusion bounds from the Higgs searches at LEP and the Tevatron. As input, it requires 
a selection of model predictions, such as Higgs masses, branching ratios, effective couplings 
and total decay widths. HiggsBounds then uses the expected and observed topological cross 
section limits from the Higgs searches to determine which points in the parameter space have 
already been excluded at 95% CL. HiggsBounds will be updated to include new results as they 
become available. 



1 Introduction 

The search for Higgs bosons is a major cornerstone of the physics programmes of past, present 
and future high energy colliders. The LEP and Tevatron experiments, in particular, have been 
able to turn the non-observation of Higgs bosons into constraints on the Higgs sector, which 
can be very useful in reducing the available parameter space of particle physics models. Such 
constraints will continue to be important far into the LHC era as they will need to be taken into 
account in the interpretation of any new physics. 

These analyses usually take one of two forms. Dedicated analyses have been carried out 
in order to constrain some of the most popular models, such as the SM^ and various bench- 
mark scenarios in the MSSJVp. In addition, model-independent limits on the cross sections of 
individual signal topologies (such as e + e~ — > h%Z — > bbZ) have been published. The former 
type of analyses include detailed knowledge of the overlap between the individual experimental 
searches, and therefore have a high sensitivity, whereas the latter can be used to test a wide 
class of models. 

There are certain issues involved with the application of these experimental constraints. The 
data is distributed over many different publications and the limits are given with a variety of 
normalisations. In the case of the Tevatron, the results are also frequently updated. Further- 
more, care must be taken when using more than one experimental analysis to ensure that the 
resulting exclusion bound has the same confidence level (CL) as each individual analysis. 

The fortran code HiggsBounds ^ has been designed to facilitate the task of comparing Higgs 

"Online version and code download available at: http://www.ippp.dur.ac.uk/HiggsBounds 
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Table 1: LEP and Tevatron analyses used by HiggsBounds. I or I' indicates an electron or a muon, and f indicates 
analyses which combine processes using SM assumptions. In this notation, mi t > m^. 
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sector predictions with existing exclusion limits, thus allowing the user to quickly and conve- 
niently check a wide variety of models against the state-of-the-art results from Higgs searches. 

2 Outline of the program 

The user provides the Higgs sector predictions of the model under consideration. For each 
neutral Higgs boson hi (i = 1, . . . , nHi ggs ) in the model, this will usually include the mass, total 
decay width, branching ratios and Higgs production cross sections: 

M ht ,T tot (hi) ,BR modcl (hi -> ...) , CT - d ;ff . (1) 

Where it exists, a SM (P) is used as the reference cross section. Variations on this input format 
are offered, as described in detail in the HiggsBounds manual^} The HiggsBounds package 
includes sample programs which demonstrate how HiggsBounds can be used in conjunction with 
the widely used MSSM Higgs sector progr ams FeynHiggs - and CPsuperH^^. 

A list of the experimental analysesPEE! currently included in HiggsBounds is given in Tabled! 
These include results from both LEP and the Tevatron and consist of tables of expected (based 
on Monte Carlo simulations with no signal) and observed 95% CL cross section limits, with 
a variety of normalisations. The list mainly consists of analyses for which model-independent 
limits were published. However, we also include some dedicated analyses carried out for the 
case of the SM. These analyses are only considered if the Higgs boson in question would appear 
sufficiently 'SM-like' to this analysis. Roughly speaking, this requires that the ratios of all 
involved couplings to the SM couplings are approximately equal^. 

For each Higgs process X (here, we treat each combination of Higgs bosons in each ex- 
perimental analysis as a separate X), HiggsBounds uses the input to calculate the quantity 
Qmodei(^Q) which, up to a normalisation factor, is the predicted cross section for X. 

The normalisation is carried out using SM predictions for Higgs boson production cross sec- 
tions and decay branching ratios from HDECAY 5 3.303, the TEV4LHC Higgs Working Group^El, 
VFBONLo'' 7 , HJET^ 1.1 and dedicated calculations of our own^U. 

In order to ensure the correct statistical interpretation of the results, it is crucial to only 
consider the experimentally observed limit for one particular X. Therefore, HiggsBounds must 
first determine Xq, which is defined as the process X with the highest statistical sensitivity 
for the model point under consideration. In order to do this, the program uses the tables of 
expected experimental limits to obtain a quantity Qoxpec corresponding to each X. The process 
with the largest value of Qmodei/Qexpec is chosen as Xq. 



Figure 1: Coverage of the LEP Higgs searches in the M/fj-tan/3 plane of the CPX scenario, where Mh x is the 
lightest neutral Higgs boson and tan (3 is the ratio of vacuum expectation values. Left: the LEP processes predicted 
to have the highest statistical sensitivity at each parameter point. Right: the parameter regions excluded at the 
95% CL (green (dark grey) = excluded, white = unexcluded, light grey = theoretically inaccessible) 
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HiggsBounds then derives a value for Qobs for this process Xq, using the appropriate table 
of experimentally observed limits. If 

Qmodcl(^o) > ^ ^) 
Qobs(^o) 

HiggsBounds concludes that this particular parameter point is excluded at 95 % CL. 

In order to use HiggsBounds, the narrow- width approximation must be valid for each Higgs 
boson described in the input. This is because the experimental exclusions bounds currently 
utilised within HiggsBounds have all been obtained under this approximation. We intend to 
include width-dependent limits into HiggsBounds in the future, where they are provided by 
experimental collaborations in a model-independent format (such as those provided in ReP). 

3 Numerical example: the CPX scenario 

We will illustrate some of the main features of HiggsBounds using an example from Ref.Q^ 
(with the modifications described in Ref.^^). The CPX scenario was one of the MSSM scenarios 
which were investigated in detail by the LEP Higgs Working Group^. It is phenomenologically 
interesting because it introduces large CP-violating phases, which induce mixing in the neutral 
Higgs sector, resulting in weaker exclusions than those obtained for the real M SSM. However, 
since this original analysis, there have been relevant theoretical advances ^ ! which can have 
a large numerical effect on the Higgs sector of the CP-violating MSSM. Therefore, HiggsBounds 
was employed to investigate the effect of these new results on the amount of CPX parameter 
space which can be excluded by current Higgs search data. 

From Fig. [1] (right), it can be seen that, although substantial regions of CPX parameter space 
can be excluded (green), there are significant regions which remain unexcluded (white), including 
an unexcluded region (the 'CPX hole') at a lightest Higgs mass M/ ll ~ 45 GeV, qualitatively 
confirming the result of the original analysis^!. In addition, the use of HiggsBounds allows 
a greater understanding of the theoretical influences on the exclusions. For example, it can 
be seen from Fig. CD (right) that the process e + e~ — > (h\Z) — > (bbZ) (red), which is usually 
the most effective at excluding areas of MSSM parameter space, has the highest statistical 
sensitivity only in regions with low tan/3 and/or high M/ ll . This is because the coupling of 
the lightest Higgs to two Z bosons is suppressed in the other regions, therefore reducing the hi 
Higgsstrahlung production cross section. It is also interesting to note that, near to the CPX 



hole, the processes with the highest statistical sensitivity all directly involve the decay of the 
second heaviest neutral Higgs /12. Therefore, it can be inferred that variations in the partial 
decay widths of the dominant decay modes (in this case, the Higgs cascade decay h% — > h\h\ and 
the decay to b-quarks — ► bb) will affect the size and position of the CPX hole, as is indeed 
the caseQ21 

In conclusion, the program HiggsBounds provides a convenient way to compare theoretical 
Higgs sector predictions with the current exclusion bounds from LEP and the Tevatron, in a 
way that maintains the statistical interpretation of the exclusion limit and gives extra insight 
into phenomenological influences on the result. 
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